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Based on the closed time path formalism, a new Feynman rule for directly calculating the retarded 
and advanced Green functions is deduced. This Feynman rule is used to calculate the two-point 
self-energy and three-point vertex correction in (f)^ theory. The generalized fluctuation-dissipation 
theorem for three-point nonlinear response function is verified. 



> 

ON 

O 

o 
o 
o 

Oh: 
I 

Oh! 

<D . 



X 



The closed time path formalism (CTPF) in real- 
time finite-temperature field theory has been used widely 
to investigate the equilibrium and non-equilibrium prop- 
erties of the thermal system recently. It is differ- 
ent to the imaginary time formalism (ITF), without ana- 
lytical continuation the physical region can be reached 
in the real time formalism. In the CTPF two kinds 
of fields are defined according to the closed time inte- 
gral contour in the generating function of Green func- 
tion. The closed time integral contour consists of the 
first branch which runs from negative infinity to posi- 
tive infinity and the second branch which runs back from 
positive infinity to negative infinity. The field located 
on the first and the second branch is defined as phys- 
ical (type-1) and ghost (type-2) field, respectively. So 
that the propagator becomes a 2 x 2 matrix with compo- 
.nents Gn, G12, G21, 6*22, their corresponding self-energy 
'are denoted as Sn, S12, S21, ^22 according to the type- 
, 1 or type-2 external leg of the Feynman diagram. It is 
'known that the self-energy with physical interpretation 
'is not Ell but the retarded function S''''* (or advanced 
• function S'^''") which is the analytical continuation of the 
self-energy obtained in the ITF . In order to calculate 
retarded (or advanced) self-energy in the real time for- 
malism, the usual way is to calculate En, E12, E21 and 
E22 first, and then use the relation S'"'^* ~ Sn + S12 
(or E"*" = Ell + E21) to get them. When one calcu- 
lates n-point {n > 2) retarded/advanced vertex for 
example, three-point fully retarded vertex Traa defined 
in Ref. we need to calculate eight components Fm, 
rii2, ri2i, F122, F211, r2i2, F221 and r222 first and then 
use the relation Traa = ^i^iii - ^112 - ^121 + ^122 + 
^211 - ^212 - ^221 + ^222) to get it, the calculation be- 
comes somewhat tedious. In this letter we deduce a novel 
Feynman rule which can be used to calculate n-point re- 
tarded/advanced Green function directly to avoid above 
tedious calculation, and then we verify the generalized 
Fluctuation-dissipation theorem for three-point nonlin- 
ear response function by calculating the vertex correc- 
tion. 

In the CTPF, for any field (f) the four components of 
the matrix propagator in the single-time representation 
are defined as 
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where Tp represents the time ordering operator along the 
closed time path, it is normal and antichronological time 
ordering on the first and second branch of the closed time 
path, respectively, ai, 02 £ {1, 2} indicate on which of the 
two branches the (f> fields are located. (• • •) stands for the 
thermal expectation value. Following |^ we define 

0a(a;) = 0i(a;)-02(a;) , ^^(a;) = |((/'i(a;)-|-(/)2(a;)) , (2) 
and the 2-point Green function in the above (r, a) basis 

Gaia2ixiyX2) = -i2"^''^^ {Tp{(j)aiixi)(l)a2{x2))) , (3) 

where ai, a2 € {a,r}, and Ur is the number of r indices 
among {ai,a2}- Inserting Eq. (||) into (||) it is easy to 
get 

Grr(a;i,a;2) = Gi2(a;i,a;2) + G2i(a;i, ^2) 

= -z([0(xi),0(x-2)]±), (4) 

Graixi,X2) = Gii{xi,X2) - Gi2{xi,X2) 

= -^0(x'^-x'2)m^l),^{^2)h), (5) 
Gar(a;i,a;2) = Gii(a;i,a;2) - G2i(a;i,a;2) 

= i0ix°-x",)mx,),<l>{x2)h), (6) 

Gaa{xi,X2)^0. (7) 

Here [4>ixi), (f){x2)]± = ^(xO^fe) ± 0(a;2)(/>(a;i). The 
double sign ± or ^ in Eq. (|j)-(|6|) corresponds to Boson 
and Fermion field for the upper and lower case, respec- 
tively. Obviously Graixi,X2) and Gar(a;i,a:2) are the 
usual retarded and advanced Green functions. 

Denote the 2x2 matrix propagator in the single-time 
representation as 

and the matrix propagator in the (r, a) representation as 

_ f Gaa Gar | _ ( ^ G"*^" \ 

' Gra Grr ' ' ' ' ^ ' 



Both matrix propagators satisfy following transformation 
relation iHl: 



G = Q^GQ, (10) 
where Q is orthogonal Keldysh transformation matrix 
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From the above transformation (|10p we see that, both 
left and right side of the propagator G associate with the 
transformation matrix Q. We can estabUsh a new Feyn- 
man rule in following way: as illustrated in Fig. |l|, the 
matrices can be absorbed into the left vertex with 
outgoing momentum line, and Q can be absorbed into 
the right vertex with incoming momentum line, so that 
we leave Gap in the (r, a) basis as a new propagator of the 
propagating line. After absorbing all Q from propagating 
line, the new bare vertex with all incoming momentum 
(as illustrated in Fig. ||(a)) can be defined as 



ll3P'-p" (p, q,---r)= Qi3aQf3'b ■ 



zgab- 



(13) 



where a, 6, • • • , c G 1,2 and /3, /?',•• • P" G a, r, a summa- 
tion over repeated indices is understood, p+q+- ■ ■+r = 
because of energy-momentum conservation. In the single- 
time representation, gii...i — — g22- -2 = 9, all other com- 
ponents vanish, here g is the coupling constant. As shown 
in Fig. 1^, we change the orientation of incoming momen- 
tum r in Fig. ^a) as outgoing line in Fig. |^(b) with 
energy- momentum conservation p + q + -- - — r = Q. Af- 
ter absorbing associated with outgoing momentum r 
and all Q associated with incoming momentum into the 
vertex, this new bare vertex in Fig. §(b) should be 



7/3/3' ■■■/3" 



{p,q, r) = QpaQ/s'b- ■ - QI 



(14) 



From E g. (pi] ) and (|T2|) we know Qj^nc ~ QI/^"- Insert it 
into Eq. (|14|) and then compare to Eq. ( p^ , we see clearly 
that the bare vertex 7^^/...^//(p, q, • • • r) with p + q+ ■ ■ ■ + 
r = is the same as bare vertex 7^/3/...^// (p, g, • • • — r) 
with p + q + - ■ ■ — r = 0. This property indicates that the 
new bare vertex defined here is independent of the orien- 
tation of the momentum. In the following we can drop 
the momentum arguments for the bare vertex. Substi- 
tuting Qal, Qa2, Qrl, Qr2 in Eq.(|ll|) iuto Eq.(|l|), we can 
express the n-point new bare vertex in the (r, a) basis as 



7q 



(-1) 



7ia,{ai,a2,---,a„)] 



(15) 



where na{ai,a2, 
among {ai, a2, • • 



• • , q;„) is the number of a indices 
, a„}. As this number is even, the bare 
vertex is zero, so that in the (r, a) basis there are a half 
of bare vertices vanish and other are only related to the 
coupling constant and independent of the thermal distri- 
bution function. This property will help us to simplify 
the calculation greatly. For example, the explicit form of 
three-point bare vertices read 



7aaa — ^arr — ^rra — 
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'Jaar — ^ara — ^rrr "~ '^raa ' ^ • 



(16) 





FIG. 1. The new bare vertex and propagator Gafi after 
absorbing and Q into the left and right vertex in our de- 
duced new Feynman rule. 





(a) 



(b) 



FIG. 2. The bare vertex in the (r, a) basis. The difference 
between (a) and (b) is that the orientation of momentum r is 
changed. 




FIG. 3. The one- loop Feynman diagram for calculating the 
retarded and advanced self-energy in cj>^ theory. 

As first application we use this new Feynman rule to 
calculate the retarded (or advanced) self-energy in (p'^ the- 
ory. As shown in Fig. 0, the self-energy in the (r, a) basis 
can be expressed as 

-tEp^,{q) = J -^i-ljaps) 

(zA_,(p))(-z7„'^',0(A,',(r)), (17) 

where r = p + q. Noticing propagator Aaa — 0, the poles 
of Ara{p)Ara{i") and Aq^ (p) A^^ (r) are both on the same 
side of the real axis in the complex plane, we have 



+00 ^ +OC ^ 

J ^Araip)Ara{r) = J ^ Aarip) Aar{r) 



= 0. 



(18) 



Then four components of the self-energy can be derived 
as 



iT,ra{q) = g 



{n{ro)Aar{p) [A,.Q(r) 



(27r)4 

+n{po)Ara{r)[Ara{p) - Aar{p)]} , 



(19) 



iT^ariq) = 5^ / ■7^;^[n'{po)Aar{r)[Ara{p) - Aar{p)] 



-iT.rr{q) 



(2^Y 

+n(ro)A™(p)[A™(r) - A^rW]} , (20) 

[2n(go) + l][S™(9)-Sar(9)], (21) 

0. (22) 



Here n(A:o) — l/[exp(/3A:o) — 1] is the Bose thermal dis- 
tribution. It is easy to show that Yiar(,q) and T^raiq) are 
the retarded and advanced self-energy which correspond 
to the retarded and advanced analytical continuation of 
the self-energy calculated in the ITF |l2). Eq.(^) corre- 
sponds to the Fluctuation-Dissipation theorem in linear 
response theory ||l^ . 

The one-loop three-point vertex correction is illus- 
trated in Fig. ^, the general expression for this vertex 
correction in the (r, a) basis can be expressed as 



iTap5{p,q,r) = 



4(-*7aa'a")[*A„'^"(/l)] 




(27r)** ^ " " VL"— a/3 

{-il0p'(3")V^l3'S"i^-2)] 
{-llss-5")[^^5'c."ih)]. 

and li+l2+h = 0. 



(23) 



FIG. 4. The Feynman diagram for calculating one-loop ver- 
tex correction in the (r, a) basis in <j!>^ theory. 

Using Aaa ~ and the Z°-integral over terms of the 

type Ara{ll)Ara{l2)Ara{h) and Aar{ll)Aar{l2)Aar{h) 

vanish because the poles of these two terms are both on 
the same side of the real axis in the complex l\ plane, we 
deduce eight components of one-loop vertex correction in 
(j)^ theory as 



dHi 



i^n(lm)[Ara{ll)-Aar{ll)\Aar{h)Ara{h) 



+n{lo2)[Ara{l2)-Aar{l2)]Ara{h)Ara{h) 
+n{hz)[Ara{h)-Aar(lz)]Aar(ll)Aar(l2)], (24) 

\^n{lQi)[Ara{liyAar{ll)]Ara{l2)Ara{k) 
+n{lo2)[Arail2)~Aar{l2)]Aar{h)Aar{k) 
+n{lo3)[Araih)-Aar{l3)]Aarih)Arail2)], (25) 
2^3 f d% 

-r_(p,g,r) = ^y ^ 

\^niloi)[Ar.a{liyAarih)]Aaril2)Aarih) 
+n{lo2)[Ara{l2)-Aar{l2)]Ara{ll)Aar{l3) 
+n{ln3)[Ara{l3)~Aar{k)]Ara{ll)Ara{l2)}, (26) 

d% 



-iTaar{p,q,r) - ^^^^j ^^^^^ 

{[Ara{ll)-Aar{h)]Ara{l2)Aar{l3)[N2-N3]Ni 

+ [ArailiyAarill)]Aaril2)Aaril3)[l~NiN2] 

+ [Ara{liyAarill)]Arail2)Arail3)[NlN3-l]}, (27) 

-lTraa{p.q,r) - J 

{[Arail2yAaril2)]Araill)Araih)[NlN2-l] 

+ [Arail2yAaril2)]Aarill)Araih)[N2N3-N2Ni] 

+ [Arail2)-Aar{l2)]Aar{ll)Aaril3)[l-N2N3]}, (28) 

dHi 



-iTaraip,q,r) - ^^^^ j ^^^^^ 

[[Arail3)-Aarih)]Aarill)Aaril2)[l~NiN3] 

+ [Ara{l3)-Aar{l3)]Ara{h)Aar{l2)[NiN3-N2N3] 

+ [Ara{l3yA,r{l3)]Ara{h)Ara{l2)[N2N3-l]Y (29) 

-.Wp,,,.)=(^/|^{A^iA^.iV3 

x[Araill)Arail2)Ara{l3)-Aarih)AaAl2)Aaril3)] 
+ {Nl-NiN2N3)[Araill)-Aarill)]Arail2)Aar(.l3) 
+ iN2-NiN2N3)[Ara{l2)-Aaril2)]Aarih)Aar{l3) 
+ iN3-NiN2N3)[Ara{l3)-Aaril3)] 

xAraih)Aaril2)}, (30) 
-irrrr{P,q,r) =0. (31) 



Here N, ^ I + 2n{hi). Eq.(||) and (|T|) indicate the 
two-point self-energy and three-point vertex with all in- 
dices being r vanish, this property is also true for n-point 
vertex in n > 3 case IKjI . 



Denote iVfeo = 1 + 2n{k°) for = p°,q°,r°. As 
p + q + r = and li + h + h — 0, we have 



Ngo[Ni - N2] = N1N2 ~ I , 
Nr.o[N2 - N3] = N2N3 - 1 . 



(32) 
(33) 
(34) 



Using above equations we can verify that three-point ver- 
tex corrections in Eq.(p4|)-(^T1) satisfy foUowing relations: 

Tara = ^p^i^rar ~ ^rra) + Nro{T*^^ — Tarr) , (35) 
Traa — ^q" (r*„. — ^rra) + ^r" (r^rr ~ ^rar) , (36) 
Taar = ^p^i^rra ~ ^rar) + Afgo{T*j,^ — ^arr) , (37) 
Taaa = r*„, + T*^^ + T*^^ + NqoNro{Ta„- + ^Irr) 
+NpoNro{Trar+Kar) 

+NpoNgo{Trra+K,,)- (38) 

These relations are just the generalized Fluctuation- 
Dissipation theorem for three-point nolinear response 
function which is in agreement with the result obtained 
in our previous work Eq.(|35|)-(|3^) together with 

Eq.(^l|) indicate that, among these vertex functions 

\^ rra^^ rar arr raa^^ ara:^ aar aaa^^ rrr} Only three 

components are independent. 

In Ref. Aurenche and Becherrawy developed a new 
Feynman rule which denoted by R, A to calculate re- 
tarded and advanced vertex function. The main differ- 
ence between their and our Feynman rule is following: In 
the {R, A) Feynman rule the matrix propagator is diag- 
onal and independent of the thermal distribution func- 
tion; the thermal distribution function is absorbed into 
the vertex, the defined new vertex depends on the relative 
orientation of the momenta and only two of the vertex 
components vanish. In our (r, a) Feynman rule the ma- 
trix propagator is non-diagonal, one component is zero 
and another component Arr is related to thermal distri- 
bution function; the defined new bare vertex is indepen- 
dent of thermal distribution function and only related to 
the coupling constant; this new bare vertex doesn't rely 
on the orientation of the momenta and one half of the 
bare vertex components vanish. So that in the practi- 
cal calculation both (i?, A) and (r, a) Feynman rule have 
their own advantage and disadvantage. We should point 
out that, in our (r, a) basis the Green function is defined 
with time-ordering and have explicit physical interpre- 
tation (see Eq.(|)-(0) and Ref. ||l§), it works in both 
coordinate and momentum space; but the (i?. A) Feyn- 
man rule, as stated by the authors themselves in Ref. 
[^, only works in the momentum space, and the {R,A) 
Green function have no explicit reference to a definite 
time-ordering. 

By suitable combination we can get the {R, A) vertex 
from the (r, a) vertex. For three-point vertex the relation 
between them can be expressed as 
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(45) 



In summary, the new Feynman rule in the (r, a) ba- 
sis is deduced from the Keldysh transformation in the 
CTPF. This Feynman rule can be used to calculate the 
retarded and advanced Green function directly. In this 
new Feynman rule, the bare vertex is only related with 
the coupling constant and independent of the orientation 
of the momenta, one half of bare vertex components van- 
ish; the new propagator depends on retarded, advanced 
propagator and the thermal distribution function. As ap- 
plication, the one-loop self-energy and three-point vertex 
correction in 0"^ theory are calculated in the (r, a) ba- 
sis, and the Fluctuation-Dissipation theorem is verified 
in the two-point and three-point cases. The difference 
between our (r, a) Feynman rule and the (i?. A) Feyn- 
man rule introduced by Aurenche and Becherrawy H is 
discussed. 
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